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JLIST  OF  TABLES 


Table 

Number 

Bearing  Size 
0.  D.  "  X  I.  D.  ' 

Speeds 
'  RPM 

Number 
of  Pads 

t 

avg 

R 

rp7, 

’  ^GR 

1 

19  X  9-1/2 

160  and  310 

6.  8  b  10 

0.  154 

;  50 

Per  Fig.  4 

2 

25  X  12-1/2 

120  and  240 

6,  8  b  10 

0.  154 

50 

Per  Fig.  4 

3 

31  X  15-1/2 

180  and  320 

6,  8  b  10 

0.  154 

50 

Per  Fig.  4 

4 

37  X  18-1/2 

180  and  320 

6,  8  b  1 0 

0.  154 

50 

Per  Fig.  4 

5 

39  X  19-1/2 

150  and  200 

6,  8  b  10 

0.  154 

50 

Per  Fig.  4 

6 

41  X  20-1/2 

100  and  200 

6,  8  b  10 

0.  154 

50 

Per  Fig.  4 

7 

45  X  22-1/2 

100  and  170 

6,  8  b  10 

0.  154 

50 

Per  Fig.  4 

8 

50  X  25 

100  and  170 

6.  8  b  1 0 

0.  154 

50 

Per  Fig.  4 

9  1 

26  X  17-1/2 

160 

8 

0.  113 

50 

Per  Fig.  47 

10 

31  X 16-1/2 

320 

8 

0.  142 

5L7 

Per  Fig.  4 

11 

35  X  18-1/2 

170 

8 

0.  133 

5L5 

Per  Fig.  4 

12 

31  X  15-1/2 

180  and  320 

6,  8  b  10 

0.  130 

50 

Per  Fig.  4 

13 

31  X  15-1/2 

180  and  320 

6,  8  b  1 0 

0.  193 

50 

Per  Fig.  4 

14 

31  X  15-1/2 

180  and  320 

6,  8  b  10 

0.  154 

53 

Per  Fig.  4 

IS 

31  X  15-1/2 

180  and  320 

6,  8  b  10 

lo.  154 

47 

Per  Fig.  4 

16 

31  X  15-1/2 

180  and  320 

8  1 

0.  154 

50 

130°F 

17 

39  X  19-1/2 

150  and  200 

8 

0.  154 

50 

130OF 

18 

19  X  9-1/2 

100 

6,  8  b  10 

0.  154 

50 

Per  Fig.  4 

19 

75  X  37-1/2 

100 

6.8  b  10 

0.  154 

50 

Per  Fig.  4 

20 

100  X  50 

100 

6,8  b  10 

0.  154 

50 

Per  Fig.  4 

21 

19  X  9-1/2 

100 

6 

0.  130 

50 

Per  Fig.  4 

21 

45  X  22-1/2 

100 

6 

0.  130 

50 

Per  Fig.  4 

21 

75  X  37-1/2 

100 

0.  130 

50 

Per  Fig.  4 

21 

100  X  50 

100 

6  ' 

0.  130 

50 

Per  Fig.  4 

22 

51-1/2  X  32 

200  and  400 

10  b  12i 

0.  117 

52.6 

Per  Fig.  4 
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LIST  OF  FIGURES 

Figure 

Number  Ceptlon 

1  Moment  Equilibrium  of  Tilting  Pad  Bearing 

2  Effect  of  Thermal  Wedge  and  Pad  Bending  on  Load  Carrying  Capacity 

3  Coordinate  Syatem 

4  Groove  Mixing  Temperature 

5  Vi acoaity- Temperature  Relation  of  2190T  Oil 

6  Meah  Notation 

7  19"  O.  D.  X  9-1/2"  1.  D.  Bearing  at  160  RPM  - 

Film  Thickneaa  Va.  Unit  Loading 

8  19"O.D.  X  9-1/2"  I.  D.  Bearing  at  310  RPM - 

Film  Thickneaa  Va.  Unit  Loading 

9  19"  O.  D.  X  9-1/2"  LD.  Bearing  at  160  RPM  - 

Maximum  Temperature  Vs.  Unit  Loading 

10  19"  O.D.  X  9-1/2"  I.  D.  Bearing  at  310  RPM - 

Maximum  Temperature  Va.  Unit  Loading 

11  19"  O.D,  X  9-1/2"  L  O.  Bearing  -  Hydrodynamic  Oil  Flow 

12  25"  O.  D.  X  12-1/2"  1.  D.  Bearing  at  120  RPM  - 

Film  Thickneaa  Va.  Unit  Loading 

13  25"  O.D.  X  12-1/2"  LD.  Bearing  at  240  RPM  - 

Film  Thickneaa  Va.  Unit  Loading 

14  25"  O.D.  X  12-1/2"  I.  D.  Bearing  at  120  RPM  - 

Maximum  Temperature  Va.  Unit  Loading 

15  25"  O.D.  X  12-1/2"  I.  D.  Bearing  at  240  RPM  - 

Maximum  Temperature  Va.  Unit  Loading 

16  25"  O.D.  X  12-1/2"  L  D.  Bearing  -  .Hydrodynamic  Oil  Flow 

17  31"  O.D.  X  15-1/2"  I.  D.  Bearing  at  180  RPM  - 

Film  Thickneaa  Va.  Unit  Loading 

18  31"  O.D.  X  15-1/2"  LD.  Bearing  at  320  RPM - 

Film  Thickneaa  Va.  Unit  Loading 

19  31"  O.D.  X  15-1/2"  LD.  Bearing  at  180  RPM  - 

Maximum  Temperature  Va.  Unit  Loading 

20  31"  O.D.  X  15-1/2"  LD.  Bearing  at  320  RPM  - 

Maximum  Temperature  Va.  Unit  Loading 

21  31"  O.D.  X  15-1/2"  1,  D.  Bearing  -  Hydrodynamic  Oil  Flow 

22  37"  O.  D.  X  18-1/2"  LD.  Bearing  at  180  RPM  - 

Film  Thickneaa  Va.  Unit  Loading 

23  37"  O.D.  X  18-1/2"  1.  D.  Bearing  at  320  RPM 

Film  Thickneaa  Va.  Unit  Loading 


Figure 

Number 


Caption 


24  37"O.D.  X  18-1/2"  I.  D.  Bearing  at  180  RPM - 

Maximum  Temperature  Vs.  Unit  Loading 

25  37"0.  D.  X  18-1/2"  I.  D.  Bearing  at  320  RPM - 

Maximum  Temperature  Vs.  Unit  Loading 

26  37"  O,  D.  X  18-1/2"  I,  D.  Bearing  -  Hydrodynamic  Oil  Flow 

27  39"  O.  D.  X  19-1/2"  1.  D.  Bearing  at  150  RPM  - 

Film  Thickness  Vs.  Unit  Loading 

28  39"  O.  D.  X  19-1/2"  I.  D.  Bearing  at  200  RPM  - 

Film  Thickness  Vs.  Unit  Loading 

29  39"0.  D.  X  19-1/2"  LD.  Bearing  at  150  RPM - 

Maximum  Temperature  Vs.  Unit  Loading 

30  39"  O.D.  X  19-1/2"  I.  D.  Bearing  at  200  RPM  - 

Maximum  Temperature  Vs.  Unit  Loading 

31  39"  O.D.  X  19-1/2"  I,  D.  Bearing  -  Hydrodynamic  Oil  Flow 

32  41"  O.D.  X  20-1/2"  I.  D.  Bearing  at  100  RPM - 

Film  Thickness  Vs.  Unit  Loading 

33  41"  O.D.  X  20-1/2"  I.  D.  Bearing  at  200  RPM  - 

Film  Thickness  Vs.  Unit  Loading 

34  41"  O.D.  X  20-1/2"  I.  D.  Bearing  at  100  RPM - 

Maximum  Temperature  Vs.  Unit  Loading 

35  41"  O.D.  X  20-1/2"  l.D.  Bearing  at  200  RPM 

Maximum  Temperature  Vs.  Unit  Loading 

36  41"  O.D.  X  20-1/2"  l.D.  Bearing  -  Hydrodynamic  Oil  Flow 

37  45"  O.D.  X  22-1/2"  I.  D.  Bearing  at  100  RPM 

Film  Thickness  Vs.  Unit  Loading 

38  45"  O.  D.  X  22-1/2"  l.D.  Bearing  atS70  RPM  - 

Film  Thickness  Vs.  Unit  Loading 

39  45"  O.D.  X  22-1/2"  l.D.  Bearing  at  100  RPM  - 

Maximum  Temperature  Vs.  Unit  Loading 

40  45"  O.D.  X  22-1/2"  l.D.  Bearing  at  170  RPM - 

Maximum  Temperature  Vs.  Unit  Loading 

41  45'  O.D.  x  22-1/2"  l.D.  Bearing  -  Hydrodynamic  Oil  Flow 

42  50"  O.  D.  X  25"  I.  D.  Bearing  at  100  RPM  - 

Film  Thickness  Vs.  Unit  Loading 

43  50"  O.  D.  x\25"  I.  D.  Bearing  at  170  RPM  - 

Film  Thickness  Vs.  Unit  Loading 

44  50"  O.D.  X  25"  I.  D,  Bearing  at  100  RPM  - 

Maximum  Temperature  Vs.  Unit  Loading 

45  50"  O.D.  X  25"  l.D.  Bearing  at  170  RPM  - 

Maximum  Temperature  Vs.  Unit  Loading 

46  50"  O.  D.  X  25"  1.  D.  Bearing  -  Hydrodynamic  Oil  Flow 
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Figure 

Number 


Ceptlon 


47 

48 

Groove  K 
26"  O.  O. 

49 

26"  O.  O. 

50 

26"  O.  O. 

51 

31"  O.  O. 

52 

31"  O.  O. 

53 

31"  0.0. 

54 

35"  O.  O. 

55 

35"  O.  O. 

56 

35"  O.O. 

57 

31"  O.O. 

58 

31"  0.0. 

59 

31"  0.0. 

bO 

31"  0.0. 

bl 

31"  0.0. 

62 

31"  0.0. 

63 

31"  0.0. 

64 

31"  0.0. 

65 

6b 

31"  0.0. 
31"  0.0. 

67 

31"  0.0. 

17-1/2"  I.  O.  Bearing  at  160  RPM  - 
Film  Tliickneas  Vs.  Unit  Loading 
17-1/2"  I.  O.  Bearing  at  160  RPM  - 
Maximum  Temperature  Vs.  Unit  Loading 

17- 1/2"  1.  D.  Bearing  at  160  RPM  - 
Hydrodynamic  Oil  Flow  Vs  Unit  Loading 

16-1/2"  I.  D.  Bearing  at  320  RPM  - 
Film  Thickness  Vs.  Unit  Loading 
16-1/2"  I.  O.  Bearing  at  320  RPM  - 
Maximum  Temperature  Vs  Unit  Loading 
16-1/2"  L  D.  Bearing  at  320  RPM  - 
Hydrodynamic  Oil  Flow  Vs.  Unit  Loading 

18- 1/2"  1.  D.  Bearing  at  170  RPM  - 
Film  Thidkneas  Vs.  Unit  Loading 

18-1/2"  l.D.  Bearing  at  170  RPM  - 
Maximum  Temperature  Vs.  Unit  Loading 
18-1/2"  1.  O.  Bearing  at  170  RPM  - 
Hydrodynamic  Oil  Flow  Vs.  Unit  Loading 
15-1/2"  LD.  Bearing  at  180  RPM, 

Film  Thickness  Vs.  Unit  Loading 
15-1/2"  1.0.  Bearing  at  320  RPM, 

Film  Thickness  Vs.  Unit  Loading 
15-1/2"  1.0.  Bearing  at  180  RPM, 


Tavg/R 

Tavg/R- 


:  0.  130 


0.  130 


T.vg/R  •  0.  130 

Maximum  Temperature  Vs.  Unit  Load* 
x  15-1/2"  LO.  Bearing  at  320  RPM,  T^^  /R  »  0.  130 
Maximum  Temperature  Vs.  Unit  Load* 

X  15-1/2"  1.  O.  Bearing  at  180  RPM,  T,,,_/R  =  0.  193 
Film  Thickness  Vs.  Unit  Loading 
X  15-1/2"  1.  O.  Bearing  at  320  RPM,  T^^  /R  •  0.  193 
Film  Thickness  Vs.  Unit  Loading 
X  15-1/2"  1.0.  Bearing  at  180  RPM,  T^^g/R^  *’5 

Maximum  Temperature  Vs.  Unit  Load&g 
X  15-1/2"  1.0.  Bearing  at  320  RPM,  •  0.  193 

Maximum  Temperature  Vs  Unit  Loading 
X  15-1/2"  I.  O.  Bearing  -  Hydrodynamic  Oil  Flow 
X  15-1/2"  1.0.  Bearing  at  180  RPM,  r  %  «  53 
Film  Thickness  Vs.  Unit  Loading 
X  15-1/2"  1.  D.  Bearing  at  320  RPM, 

Film  Thickness  Vs.  Unit  Loading 


rp%  - 


53 
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Figure 

Number 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 

81 

82 

83 

84 

85 

8b 

87 

88 

89 

90 


Caption 


31"  O.  D  X  15-1/2"  I.  D.  Bearing  at  180  RPM,  rp%  *  53 
Maximum  Temperature  Vs.  Unit  Loading 
31"O.D.  X  15-1/2"  I.  D.  Bearing  at  320  RPM,  rp%  •  53 
Maximum  Temperature  Vs.  Unit  Loading 
31"0.  D.  X  15-1/2"  1.  D.  Bearing  at  180  RPM,  rp%  «  47 
Film  Thickness  Vs.  Unit  Loading 
31"0.  D.  X  15-1/2"  I.  D.  Bearing  at  320  RPM,  rp%-«  47 
Film  Thickness  Vs.  Unit  Loading 
31"0.  D.  X  15-1/2"  1.  D.  Bearing  at  180  RPM,  rp%  -  47 
Maximum  Temperature  Vs  Unit  Loading 


31"  O.  D.  X 


15-1/2"!.  D.  Bearing  at  320  RPM,  rp%  »  47 
Maximum  Temperature  Vs.  Unit  Loading 
15-1/2"  1,  D.  Bearing  -  Hydrodynamic  Oil  Flow 
15-1/2"  I.  D.  Bearing,  Tq^^  =  130°F 
Minimum  Film  Thickness  Vs.  Unit  Loading 
31"O.D.  X  15-1/2"  I.  D.  Bearing,  Tf-,^  i  130OF 

Max-mum  Temperature  Vs.  Unit  Loading 


31 

31 


O.  D. 
O.  D. 


1.  D.  Bearing.  Tqj^  *  UO^F 


31"  O.D.  X  15-1/ 

Hydr.  ynamic  Oil  Flow  Vs.  Unit  Loading 
39"  O.D.  X  19-1/2"  I.D,  Bearing.  T^p  .  130OF 

Minimum  Film  ThlcknessVs.  Unit  Loading 
39"  O.D.  X  19-1/2"  I.D.  Bearing,  Tq^  :  130®F 

Maximum  Temperature  Vs.  Unit  Loading 
39"  O.D.  X  19-1/2"  I.D.  Bearing,  Tqj^  =  130°F 

Hydrodynamic  Oil  Flow  Vs.  Unit  Loading 
Bearing  Size  Vs.  Unit  Loading 
51-1/2"  O.  D.  X  32"  I.  D,  Bearing  at  200  RPM  - 

•linimurn  Film  Thickness  Vs.  Unit  Loading 
O.  O  '  32"  1  D.  Bearing  at  200  RPM  - 

Maximum  Temperature  Vs  Unit  Loading 
O.D.  X  32'  ID  Bearing  at  400  RPM- 

Minimum  Film  Thickness  Vs.  Unit  Loading 
O.  U.  X  i,  D.  ...  V/  KPM  - 

Maximum  Temperature  Vs.  Unit  Loading 
51-1/2"  O.D.  X  32"  I.D.  Bearing  -  Hvdrodynamic  Oil  Flow 
Chart  of  Maximum  and  Average  Temperatures 
Chart  of  Load  Carrying  Capacity,  6  pads,  h.. 

Chart  of  Load  Carrying  Capacity.  6  pads. 

Chart  of  Load  Carrying  Capacity,  6  pads,  h' 


51-1/2' 

51-1/2' 
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BACKGROUND  AND  SCOPE 


In  1958  the  General  Engineering  Laboratory  made  a  study  of  propeller  shaft 
thrust  bearing  operation  and  reported  its  findings  in  Reference  1.  Following  this 
study  a  comprehensive  analytical  and  experimental  program  was  undertaken,  for 
the  purpose  of  extending  present  understanding  of  these  bearings  and  in  order  to 
provide  a  body  of  design  information  for  use  in  bearing  design  and  selection.  This 
program,  like  the  preceding  introductory  study,  is  being  performed  under  a 
contract  awarded  by  the  Bureau  of  Ships  to  General  Electric  Company's  Medium 
Steam  Turbine,  Generator  and  Gear  Department. 

The  program  is  divided  into  three  phases  as  follows: 

Phase  I :  Investigate  analytically  the  performance  of  propeller  shaft  thrust 

bearings  using  the  existing  Reynolds- Energy  Method  of  solution 
to  provide  data  necessary  in  design  and  selection  of  these  bearings. 

Phase  II:  Extend  existing  analytical  techniques  for  propeller  shaft  thrust 

bearings  by  including  a  numerical  method  of  solutioi  of  the  Elasticity 
Equation.  Review,  and  where  necessary,  modify  the  design  data 
obtained  in  Phase  1,  so  as  to  include  the  effects  of  pad  distortion 
caused  by  pressure  distribution  and  thermal  gradients. 

Phase  III:  Instrument  thrust  bearings  on  two  U,  S,  Navy  ships  and  obtain 

experimental  data  on  the  performance  of  these  bearings.  This  data 
is  to  be  obtained  in  tests  carried  out  at  the  time  of  scheduled  sea 
trials.  The  thrust  bearing  performance  measurements  obtained  in 
these  sea  trials  is  to  be  used  for  correlation  with  the  design  data 
obtained  analytically  in  Phases  I  and  II. 

A  fourth  phase  which  included  the  building  of  a  thrust  bearing  test  stand  was  con¬ 
templated  but  was  not  included  in  the  present  program,  since  the  findings  of  the 
program  could  be  used  to  determine  the  features  of  the  stand. 

The  following  is  our  final  report  on  Phase  I. 
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I  INTRODUCTION 

Virtually  all  ships  in  this  country,  both  merchant  and  combat,  use  Kingsbury 
type  tilting  pad  bearings  to  transmit  the  propeller  thrust  to  the  hull  of  the  vessel. 
The  geometry  of  these  bearings  and  the  principles  on  which  they  operate  are  well 
known  and  are  described  in  most  text*  on  lubrication  as  well  as  in  the  catalogues 
of  Kingsbury  Machine  Works. 

Generally,  the  bearing  pads  are  centrally  pivoted,  i.  e. ,  the  spherical 
pivot  back  of  each  pad  is  located  mid-way  between  the  leading  and  the  trailing  edges. 
Central  pivot  location  is  required  for  reversibility,  i.  e. ,  for  operation  under  either 
direction  of  shaft  rotation. 

Conventional  bearing  calculations  in  which  temperature  variations  in  the 
oil  film  are  neglected  and  in  which  a  converging  wedge  is  formed  by  the  tilting 
of  a  flat  pad  fail  to  predict  the  load  carrying  capacity  of  centrally  pivoted  bearings. 
The  reason  for  this  is  illustrated  in  Figure  1  (a)  which  shows  the  hydrodynamic 
pressure  profile  that  is  generated  between  flat  surfaces  separated  by  fluid  film  that 
converges  slightly  in  the  direction  of  motion.  Under  these  conditions,  calculations 
show  that  the  resultant  of  the  hydrodynamic  pressures  lies  downstream  of  the  radial 
centerline  of  the  pad.  Since  the  reaction  to  these  pressures  must  pass  through 
the  pivot,  a  moment  exists  which  tends  to  eliminate  the  convergence  and  hence 
load  carrying  capacity.  However,  when  the  temperature  variations  in  the  oil  film 
and  the  deformation  of  the  pad  under  load  are  considered  in  the  analysis,  the 
somewhat  paradoxical  result  obtained  above  is  eliminated.  Calculations  then  show 
that  there  is  a  value  of  pad  inclination  (generally  other  than  zero)  for  which  the 
resultant  of  the  hydrodynamic  pressures  passes  through  the  central  pivot  as  shown 
in  Figure  1  (b). 

As  the  oil  flows  through  the  bearing  gap,  its  temperature  rises  due  to  the 
shearing  of  the  film.  This  rise  in  temperature  produces  a  "thermal  wedge"  action 
which  accounts  for  part  of  the  load  carrying  capacity  of  the  bearing  (provided  that 
the  viscosity  and  mass  density  of  the  lubricant  decrease  with  temperature, rise^ which 
is  the  case  for  all  known  oils).  Early  in  the  program,  calculations  were  made  to 
determine  the  magnitude  of  the  thermal  wedge  effect  in  a  centrally  pivoted  finite  pad. 
A  sector  shaped  pad  of  a  31"  8- shoe  bearing  was  analyzed,  first  with  the  pivot  in 
optimum  position  and  then  with  the  pivot  centrally  located.  In  both  cases,  the  pad 
was  assumed  to  remain  flat  and  the  other  operating  conditions  were: 

Speed  -  320  RPM 

Minimum  Film  Thickness  -  0.  001" 

Oil  -  2190T 

Oil  Temperature  at  Pad  Inlet  - 


noV 
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The  results  are  shown  in  Figures  2  (a)  and  2  (b>.  It  is  seen  that  the  thermal  wedge 
effect  allows  the  flat  centrally  pivoted  pad  to  carry  approximately  56%  of  the  load 
carried  by  toe  pad  with  optimum  pivot.  At  the  same  time,  the  maximum  temperature 
reached  with  the  flat  centrally  pivoted  pads  is  25°F  higher  than  that  reached  in  the 
pad  with  optimum  pivot.  Experience,  however,  suggests  that  the  difference  in 
performance  between  central  and  optimum  pivot  locations  is  not  so  severe.  It  was. 
toerefore,  decided  at  that  time  that,  in  order  to  make  a  more  realistic  analysis  in 
Phase  I,  it  should  be  extended  to  incorporate  a  simplified  elasticity  approach  which 
allows  pad  deformation  to  be  approximated  and  included  in  the  calculation.  Figure 
2  (c)  shows  toe  load  and  maximum  temperature  rf  the  31”  S-shoe  centrally 
pivoted  bearing  pad  under  toe  same  operating  conditions  but  with  pad  deformation 
included.  Comparison  of  Figures  2  (a)  and  2  (c)  now  show  that  the  centrally 
pivoted  pad  is  capable  of  carrying  approximately  92%  d  the  load  carried  by  the  Hat  pad 
with  optimum  pivot.  Its  maximum  temperature  is  b^F  higher  than  that  of  the  flat 
pad  with  optimum  pivot.  These  results  are  in  better  agreement  with  experimental 
evidence  and  toe  method  of  analysis  which  includes  a  simplified  elasticity  solution 
has  been  used  in  all  succeeding  calculations. 

{To  the  extent  that  pad  deformation  was  included  in  the  Phase  1  calculations 
toe  results  presented  in  this  report  have  anticipated  those  to  be  obtained  in  Phase 
II.  In  the  latter  phase,  the  Elasticity  Equation  is  more  rigorously  sdved  and 
includes,  in  addition,  thermal  deformation  of  the  pad.  However,  it  requires  a 
considerably  more  elaborate  digital  computer  program.  It  may  be  expected  that 
comparison  between  the  two  sets  of  results  will  suggest  modifications  of  the  Phase  I 
approach  to  yield  a  simple  yet  suffic.ently  rigorous  method  c£  solution.  ) 

The  conflict  between  the  isothermal,  flat  pad  method  of  sdution  and 
experience  with  entrally  pivoted  pads  has  been  realized,  since  the  time  that  Albert 
Kingsbury  accomplished  his  pioneering  work  on  slider  bearing  performance  (Ref.  2). 
More  recently,  interest  in  the  effects  of  thermal  wedge  and  of  pad  deformation  has 
resulted  in  analytical  studies  of  infinitely  wide  bearings,  some  of  which  are 
reported  in  References  3,  4  and  5.  For  the  case  of  the  finite  bearings,  the 
importance  of  including  the  effects  of  temperature  variations  in  the  oil  film  has 
been  studied  by  one  of  the  authors  of  this  report  and  it  is  explained  in  Reference  6. 

To  toe  best  knowledge  of  the  authors,  toe  present  report  is  the  first  published  study 
of  finite  centrally  pivded  pad  bearings  in  which  the  effects  of  radial  and  tangential 
inclinations,  temperature  variatiois  in  the  oil  film,  and  pad  deformation  are  all 
considered.  The  results  obtained  have  shown  good  agreement  with  experience. 

They  have  indicated  that  pad  deformations  are  of  the  order  of  the  minimum  film 
thicknesses  and  they  have  explained  such  test  results  as: 

1.  bearing  failures  caused  by  high  pad  temperatures. 
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2.  occurrence  of  bearing  failures  in  the  vicinity  of  the  pivot. 

3.  insensitivity  of  trailing  edge  film  thicknesses  at  high  loads. 

In  order  to  make  the  Phase  I  study  as  complete  as  possible,  approximatelv 
70%  more  cases  were  analyzed  than  were  called  for  in  the  coitract  for  this  phase. 

In  all,  262  operating  points  were  calculated.  At  each  operating  point,  the  values 
of  radial  and  tangential  pad  inclination  which  satisfied  equilibrium  of  moments  were 
obtained  using  a  trial  and  error  procedure.  This  procedure  required  an  average 
of  5  solutions  of  the  Reynolds  and  Energy  Equations  for  each  operating  point,  so  th.at 
the  total  number  of  solutions  exceeded  1300. 

The  studies  were  conducted  as  follows: 

1.  Eight  standard  bearings  were  analyzed  which  scanned  the  range  of  present 
day  propeller  shaft  bearing  sizes  (19"  O.  D.  to  50"  O.  D.  )  and  propeller 
speeds  (100  R.  P.  M.  to  320  R.  P.  M.  ).  Each  bearing  was  analyzed  at  two 
speeds  and  with  6,  8  and  10  pad  geometries.  Calculations  were  made  at 
three  values  of  minimum  film  thicknesses,  at  each  speed  and  geometry. 

These  calculations  have  yielded  the  value  and  location  of  the  minimum  film 
thickness,  the  temperature  and  pressure  distributions,  the  oil  flow  and 
horsepower  loss  as  functions  of  bearing  size,  number  of  pads,  unit  load 
and  speed.  In  particular,  they  have  shown  the  optimum  number  of  pads 
as  a  function  of  bearing  size,  unit  load  and  speed. 

2.  The  ahead  bearing  of  D0933  (U.  S.  S.  Barry)  was  analyzed  at  the  full  speed 
ahead  conditions.  Its  astern  bearing  was  similarly  analyzed  at  full  speed 
astern  condition.  (The  U.S.  S,  Barry  was  earlier  selected  by  the  Bureau  ol 
Ships  to  be  the  first  ship  (or  sea  trial  thrust  bearing  tests  under  Phase  III 

of  the  program.  The  thrust  bearings  of  the  starboard  shaft  of  this  ship  were 
instrumented  and  the  tests  at  sea  have  just  been  completed. ) 

3.  The  ahead  bearing  of  OLl  (U.S.  S.  Norfolk)  was  analyzed  under  full  speed 
ahead  operation.  This  bearing  was  selected  for  analysis  because  of  the 
past  history  of  several  successive  failures. 

4.  A  31"  O.  D.  X  15-1/2"  I.  D.  bearing  was  extensively  analyzed  in  order  to 
investigate  the  effects  of  pad  thickness  and  radial  pivot  location  on  bearing 
performance. 

5.  A31"O.D.  X  15-1/2"  I.O.  and  the  39"  O.  D.  x  19-1/2"  I. D.  bearing  were 
further  analyzed  to  determine  the  effect  of  pad  inlet  oil  temperature  on 
bearing  performance.  In  particular  the  effect  of  pad  inlet  temperature  on 
load  carrying  capacity  and  maximtim  temperature  were  investigated. 
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Additional  bearing  sizes  ranging  up  to  100”  O.  D.  were  analyzed  to  investigate 
the  relationship  between  bearing  size,  pad  thickness  and  unit  loading. 


^  O.  D.  X  32"  1.  D.  bearing  with  10  and  12  shoe  geometries  was 

analyzed  at  200  R.  P.  M.  and  at  400  R.  P.  M.  These  analyses  were  made 
under  separate  contract  with  General  Electric  Company's  Medium  Steam 
Turbine,  Generator  and  Gear  Department  who  authorized  their  inclusion 
ir  kne  present  report.  They  are  of  interest  because  the  upper  speed  is  quite 
high  and  the  results  illustrate  some  of  the  thrust  bearing  operating 
characteristics  that  wrill  be  encountered  as  propeller  speeds  are  raised. 
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n  ANALYSIS  AND  METHOD  OF  CALCULATION 


The  important  coneiderationt  in  thruat  bearing  analyaie  are: 

1.  Preasure  diatributioa  and  hence  load  carrying  capacity 

2.  Temperature  diatribution 

3.  Location  of  the  center  of  preaaure 

4.  Oil  flow 

5.  Horaepower  loaa 

all  aa  functiona  of  the  bearing  geometry,  film  ahape  and  apeed. 

In  thia  section,  the  equationa  from  which  theae  quantitiea  can  be  calculated 
are  given.  The  film  ahape  which  includea  the  effecta  of  pad  inclinationa  and 
deformation  is  discussed,  aa  are  the  groove  mixing  temperature  and  the  viscosity- 
temperature  relati^in.  Before  proceeding  to  these,  however,  it  is  necessary  to 
point  out  here  the  principal  limitations  of  the  analysis. 

1.  The  analysis  applies  only  to  steady  state  conditions.  It  does  not  supply 
any  information  on  the  transient  conditions  that  occur  during  start  up  and  shut 
down.  It  also  does  not  apply  to  dynamic  load  conditions  (such  as  crash  ahead 
and  crash  astern)  when  relative  axial  motion  between  the  runner  and  the  bearing 
introduces  squeese  film  effects. 

2.  Laminar  conditions  prevail  in  the  oil  film.  Actually  the  Reynolds  Number 
in  present  day  propeller  shaft  thrust  bearings  is  small  enough  for  this  condition  to 
be  satisfied  under  steady  state  operation.  This  is  illustrated  by  the  following 
calculation  for  an  extreme  case; 


U 


D  ■  50" 

N  ■  320  RPM 
h  :  0.  002" 

■  1  X  10*^  lb. -sec  /in.  2 

=  0.803  X  10*^  lb. -sec.2/in.< 
max  »  840  in. /sec. 


Reynolds  Number  (Maximum) 


135 


3.  Oil  inertia  effects  are  negligible.  At  the  relatively  low  surface  speeds  of 
propeller  shaft  thrust  bearings,  this  assumption  too  is  quite  valid. 

4.  The  fluid  is  incompressible. 

5.  Variation  of  the  specific  heat  of  the  oil  with  temperature  are  neglected. 
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A.  Reynolds  and  Energy  Equations  and  Their  Boundary  Conditions 

The  Reynolds  Equation  describes  the  hydrodynamic  pressures  generated  in 
the  oil  film  of  a  bearing.  These  pressures  separate  the  bearing  and  runner  surfaces 
when  there  is  a  relative  motion  between  them.  For  a  finite  pad,  the  Reynolds 
Equation  in  polar  coordinates  is  (Ref.  7); 


9  /rh3  _i2\  _^/h£  3p] 
ir  \  ^  Sr  )*  \jjir  30/ 


6  (Or 


ah 

ae 


(1) 


The  boundary  conditions  that  are  needed  for  the  solutim  of  this  equation  arise 
from  the  fact  that  the  pressure  falls  to  zero  at  the  pad  perimeter. 


With  the  coordinate  system  shown  in  Figure  3,  the  bcandary  conditions  are 

then: 


P  =  P  =  P  =  P  =0  (2) 

(r.  0)  (R-L,  0)  (r,6^)  (R,  0) 

Because  the  oil  film  may  break  down  in  diverging  regions  in  the  bearing,  it  is 
necessary  to  impose  an  additional  condition  which  states  that  the  pad  pressures 
never  fall  below  atmospheric. 


In  order  to  include  in  the  analysis  the  effects  of  temperature  (and  hence 
viscosity)  variations  in  the  oil  film,  the  Energy  Equation  has  to  be  solved  together 
with  the  Reynolds  Equation.  The  Energy  Equation  can  be  written  (Ref.  7); 
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12^ 


3  p  ^ 
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-  Cppg  J 

/r  CJ  h  h^ 

3  P'1  9  T 

h3 
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12yx 

dr  9r 
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xn  Equation  (3)  it  is  assumed  that  the  oil  flow  through  the  clearance  space  is 
adiabatic.  All  the  heat  generated  within  the  fluid  due  to  fluid  shear  is  considered 
to  be  carried  away  by  the  mass  transfer  of  the  fluid  and  no  heat  is  gained,  or  lost 
through  the  bearing  surfaces.  This  is  a  comparatively  good  assumption,  for  the 
heat  transfer  coefficient  at  the  fluid  boundaries  is  very  small.  (Reference  8) 


The  boundary  conditions  used  for  the  solution  of  Equation  3  are  that 


a)  the  pad  inlet  oil  is  at  the  groove  temperature  and 

b)  the  radial  temperature  gradient  is  zero  along  the  inner  and 
outer  circumferences  to  the  pad  because  of  the  cooling  effect  of  the 
surrounding  oil. 


(3) 
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Thus: 


(r,  o) 


GR 


DT 
a  r 


(R-L  6) 


II 

hr 


(R.  0) 


(4) 


With  the  introduction  of  the  proper  film  shape,  the  solutions  of  Equations  (1)  and 
(3)  yield  the  pressure  and  temperature  distributions  on  the  bearing  pads. 

B.  Load  Carrying  Capacity 

The  total  reaction  of  each  bearing  pad  and  hence  the  load  it  carries  is 
given  by  the  integral  of  the  hydrodynamic  pressures  over  the  pad  area.  Thus; 

R.  0  'P 

W  a  A  y  p  r  dr  de  (5) 

^R-L  ^  0 


C.  Oil  Flow 


Oil  is  introduced  into  each  pad  through  the  clearance  space  at  its  leading 
edge,  by  the  motion  of  the  runner.  Part  of  this  oil  leaves  the  clearance  space 
in  the  same  manner  from  the  trailing  edge.  The  remaining  part  of  the  oil  is  forced 
out  from  all  edges  by  the  pressure  gradients  that  are  built  up  over  the  bearing 
surface.  Referring  to  Figure  3,  the  oil  flow  (in  O.  P.  M.  )  through  the  four  edges 
is : 


Flow  into  the  pad: 

R 

j 

k-i  '  -  /e:o 

Flow  out  of  the  pad: 


231  Q,  .  /  [Wrh 

60  *  ■  L.  .  \  2 
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dr 
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(6) 
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r  :  R-L 
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R  R 


0  r=R 


It  is  seen  that  the  flow  through  the  leading  and  trailing  edges  it  made  of  two 
components.  The  first  of  these  is  independent  of  the  pressure  gradients  and  it 
is  referred  to  as  the  "Shear  Flow".  The  second  component  depends  cn  the  pressure 
gradients  and  it  is  referred  to  as  the  "Pressure  Gradient  Blow". 

Since  there  is  no  relative  radial  motion  between  the  runner  and  the  bearing 
pads,  there  is  only  "Pressure  Gradient  Flow"  out  of  the  inner  and  outer  circum¬ 
ferential  boundaries  of  the  bearing  pads. 

D.  Horsepower  Loss 

It  is  assumed  that  all  of  the  heat  generated  in  the  oil  film  goes  into  temperature 
rise  of  the  oil.  Thus,  the  horsepower  loss  can  be  computed  frcm  the  oil  flow 
through  the  bearing  pads  and  its  temperature  rise.  Thus: 


i  the  total  horsepower  loss  per  pad  is 
HP  =  HP2  +  HP3  HP4 


(9) 
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E.  Center  of  Pressure 


The  point  on  the  pad  surface  through  which  the  resultant  of  the  hydrodynamic 
pressures  acts  is  called  the  center  of  pressure.  Its  coordinates  are  given  by; 


r 

CP 


e 

cp 


(K  , 

2 

®T 

2' 

)l/2 

j  y  y  p  r^  cos  0  dr  de 

♦ 

y  y  p  r^  sin  0  dr  de 

(Ir-l  0 

.R-L  0 

w 


r  R  ^ 

I  j  P  r^  sin  0  dr  de 
tan-  _ 

/  /  p  r^  cos  0  dr  d0 

jk-L  0 


(10) 


F,  Film  Shape 

Under  the  hydrodynamic  pressures  and  the  pivot  reaction,  each  pad 
bends  so  that  the  bearing  surface  becomes  slightly  convex,  as  shown  m  Figure 
1  (b).  The  shape  that  the  pad  assumes  under  load  can  be  calculated  from  a 
solution  of  the  Elasticity  Equation  and  this  is  done  in  Phase  II.  For  the  present, 
however,  it  is  assumed  that  the  bearing  surface  becomes  very  slightly  spherical. 

In  accordance  with  plate  theory,  the  bending  deflections  are  taken  to  be  prop'.rtio ml 
to  load  and  inversely  proportional  to  the  pad  thickness  cubed.  Since  the  pads  are 
ball  seated,  they  also  tilt  in  both  radial  and  tangential  directions,  till  moment 
equilibrium  is  satisfied. 

The  film  shape  is  then  (see  Appendix): 


h  »  ha  mg  l^r^ 
1 


®T  ®  ' 

sin  (fij  -  -j”  )  -  r  sin  (0  -  ) 


0  r. 

T  'I' 

r^  cos  '"2"^  ” 


r2  .  r^2  -  2  r  rp  cos  (0  -  0p)+2  r^  rp  cos  (0^^  -  0p) 


M  i 
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For  cases  where  loads  are  light  and  the  bending  deflections  are  small,  it  is 
convenient  to  use  as  reference,  the  point  at  the  inside  radius  and  trailing  edge 
of  the  pad.  Equation  (11)  then  becomes  (for  a  centrally  pivoted  pad); 


h  =  hj+m.  (R-L)  sin  ~2~  -  f  ‘in 


6  "I  r  0  *  ’ 

(0  -  -  m^.  (R-Li)  qn*  ^  -  f  «0I  (tl  -  ^  ) 

J  L  J 


2Rc 


-  (R  -L)  +  2(R-L)  tp  cos  -  2  r  r  cos  (6  -  -j-  ] 


(12) 


For  cases  where  the  loads  are  large  and  the  bending  deflections  are  of  the  same 
order  as  the  minimum  film  th.ckness,  the  point  of  minimum  film  thickness  may 
fall  within  the  pad  boundary.  It  is  then  more  convenient  to  use  this  point  as  reference. 
Its  coordinates  can  be  obtained  by  differentiating  Equation  (11)  and  setting; 


o  h 
8  r 


=  0  and  —  •  0 
d  d 

The  coordinates  of  the  point  of  minimum  film  thickness  are  then  found  to  be: 


(13) 


Substituting  Equation  (13)  into  Equation  (11),  the  film  thickness  profile  becomes 
(for  a  centrally  pivoted  pad): 


m2 

mg 


+  r 


sin  (6  ~)*(^^  -  m,)  cos  (0  -  -^)  -  ^ 


In  Equations  11  through  14  above,  R^  is  the  radius  of  curvature  of  the  bent  pad 
In  the  present  analysis,  it  is  calculated  from  the  load  and  pad  thickness  (see 
Appendix)  using  the  relation: 


(14) 


0.  75  X  10-8 

t  avg 


where  W  is  the  load  per  pad 

^avg  mean  thickness  of  the  pad 


(15) 


Oil  Groove  Temperature 


Thfi  temperature  at  which  the  oil  enters  the  clearance  space  between  the 
runner  and  the  pad*  hae  an  imporUnt  effect  on  the  load  carrying  capacity  of  the 
bearing.  It  is  introduced  in  the  analysis  as  one  of  the  boundary  conditions  of 
Equation  (3). 

In  general,  the  temperature  of  the  oil  in  the  feed  grooves  between  the  pads 
is  several  degrees  higher  than  at  the  housing  inlet  ports.  This  difference  is 
largely  due  to  the  mixing  in  each  groove  with  hot  oil  discharged  from  the  trailing 
edge  of  the  downstream  pad.  It  is,  therefore,  significantly  affected  by  such 
factors  as; 

a)  quantity  of  oil  admitted  to  the  housing  (this  is  generally  several 

times  the  amovint  that  flows  through  the  clearance  spaces.  ) 

b)  extent  of  the  grooves 

c)  pad  discharge  temperature 

In  the  Phase  I  calculations,  the  pad  groove  temperatures  are  obtained  from 
the  experimental  data  of  several  investigators.  The  experimental  points  are 
plotted  in  Figure  4  and  a  represenUtive  curve  is  drawn  through  them.  This  curve 
shows  the  feed  groove  temperature  at  a  function  of  the  unit  load  carried  by  the 
bearing,  vidten  the  oil  temperature  at  housing  inlet  is  115  F. 

Figure  4  it,  of  course,  an  average  curve.  In  the  experimental  work  on 
udiich  it  is  based,  the  oil  flow  through  the  bearing  housing  was  four  to  five  times 
the  clearance  flow  and  the  total  area  of  the  grooves  was  15%  of  the  effective 
runner  area. 

Different  values  of  these  quantities  or  the  location  of  major  heat  sources  or 
sinks  near  the  bearing  housing  would  be  expected  to  affect  the  groove  temperature. 

H,  Viscosity- Temperature  Relation 

The  viscosity-temperature  relation  of  the  lubricant  is  required  in  the 
simultaneous  solution  of  Equations  (l)and  (3).  In  all  the  Phase  I  calculations, 
the  lubricant  properties  used  were  those  of  2190T  oil. 

The  absolute  viscosity  versus  temperature  plot  for  2190T  oil  is  shown  in 
Figure  5, 
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I.  Numerical  Solution  of  the  Reynolds  and  Energy  Equation! 

A  finite  difference  procedure  Mrae  used  to  solve  the  Reynolds  and  Energy 
Equations.  These,  however,  were  first  put  In  dimensionless  form  (Equations  A- 2 
and  A- 5  of  the  Appendix)  In  order  to  facilitate  comparison  between  geometrically 
similar  bearings. 


The  finite  difference  form  of  the  dimensionless  Rsynolds  and  Ensrgy 
Equations  are  given  by  Equations  A-4  and  A-7  of  the  Appendix.  Thess  are  two 
sets  of  algebraic  equations  that  can  be  solved  on  a  digital  computer  using  an 
Iterative  procedure.  Their  solutions  yield  the  pressure  and  temperature  profiles 
over  the  pad  surface. 


Figure  6  is  a  typical  thrust  bearing  sector  pad,  divided  Into  a  mash  of 
m  X  n  sections.  Referring  to  Equation.  A-4  and  Figure  6,  It  Is  seen  that  the 
dimensionless  pressure  n,J  »ay  point  Is  expressed  la  terms  of  the  corresponding 
dimensionless  pressures,  viscosities  and  film  thicknesses.  The  boundary 
condition  states  that  the  pressure  Is  aero  around  the  periphery  of  the  sector. 

In  order  to  meet  this  condition,  the  pressures  at  fictitious  Image  points  outside 
the  boundary  are  set  equal  in  magnitude  but  opposite  in  sign  to  the  pressures  at 
the  corresponding  points  Inside  the  boundary.  By  employing  a  process  of  Iteration 
the  m  X  n  equatlon^represented  by  Equation  A-4  are  solved  on  the  computer 
and  the  pressures  pi,  j  are  determined  at  each  mesh  caaterpolnt.  The  process 
of  iteration  is  continued  until  .the  difference  between  successive  values  of  the 
sum  of  the  preesures'converges  to  within  a  prescribed  error.  In  this  analysis, 
the  error  is  specified  to  be  less  than  0. 1%,  1.  e. 


Error 


r  i: 

j**i  1=1  L 

m  n 


'Pi.jJk 


KiK-i 


<  0  001 


j=l  i«l 


(16) 


The  load  carrying  capacity  of  a  bearing  is  greatly  influenced  by  the  oil 
viscosity.  The  tmnperature  (hence  the  viscosities)  at  each  mesh  point  are 
obtained  from  the  .solution  of  Equation  A- 7  .  The  boundary  conditions  for  this 
equation  are  introduced  by  setting:  a)  the  temperature,  along  tiie  Inlet  edge 
equal  to  the  groove  temperature  and  b)  the  temperatures  at  Retittous  image 
points  outside  the  inner  and  outer  circumferential  boundaries  equal  in  magnitude 
and  sign  to  those  at  the  corresponding  points  inside  die  boundaries. 

The  steps  for  the  simultaneous  solution  of  the  Reynolds  and  Energy  Equatior  -^ 
are  then  performed  in  the  following  msumer: 
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(1 )  The  value  of  the  film  thickness  at  every  point  is  determined. 

(2)  pj  .  is  assumed  equal  to  zero  and  the  known  value  of  inlet  temperature 
is’assigned  to  Ti  > 

(3)  The  values  of  j  are  then  determined  at  every  point  from  equation  A- 7 

(4)  The  values  of/^j  j  ar_e  calculated  from  values  of  Tj^  j 

(5)  Having  the  values  of/'^  ,,  h.  .  and  p^  the  first  approximation  of  the 
pressure  field  is  determinea'nrom  equation  A-4  and  improved  several 
times  by  iteration. 

(6)  The  value  of  the  pressure  field  thus  obtained  is  used  to  recalculate  the 
temperature  distribution  from  which  a  new  set  of/*^^  j  values  is  determined. 

(7)  A  second  approximation  of  the  pressure  field  is  now  obtained.  This  cycle 
of  pressure  and  temperature  iterations  is  continued  until  the  error,  which 
is  the  difference  between  successive  values  of  the  pressure  field,  falls 
within  the  limit  prescribed  in  equation  16. 

(8)  The  final  value  of  the  pressure  field  is  then  used  to  compute  the  final 
value  of  the  temperature  field. 

(9)  The  toUl  pad  load,  oil  flow,  horsepower  loss  and  the  coordinates  of  the 
center  of  pressure  are  calculated  by  means  of  Equations  A- 8,  A- 10  , 

A- 15  and  A- 16. 


j_  Trial  and  Error  Procedure  for  Pivoted  Pad  Bearings 

At  each  operating  point,  the  pad  deformation  has  to  be  related  to  the  pad 
load  in  accordance  with  Equation  15.  The  film  shape  which  depends  on  this 
deformation  and  on  the  inclinations  of  the  pad  has  to  be  such  that  the  resulting 
center  of  pressure  passes  through  the  pivot.  Finally,  the  groove  temperature 
used  has  to  be  related  to  the  unit  loading  in  accordance  with  Figure  4,  In  order 
to  meet  these  requirements,  the  following  trial  and  error  procedure  was  used: 

(1 )  For  the  bearing  geometry  being  studied,  select  a  value  of  minimum  film 
thickness. 

(2)  Estimate  the  corresponding  \mit  load  and  hence  the  groove  temperature 
(Tqo)  and  the  bending  coefficient  (K  «  l/2Rg). 

(3)  Select  values  of  radial  and  Ungentlal  inclinations  (m^  and  respectively). 

(4)  Introduce  the  above  as  input  data  and  obtain  the  corresponding  computer 
solution. 

(5)  From  the  computer  output  data  determine  the  coordinates  of  the  center 

of  pressure  and  the  actual  unit  load  (and  hence  the  actual  groove  temperature 
and  bending  coefficient).  Check  whether  these  agree  with  the  estimated 
ones  within  the  following  error  limits: 

a) 

b) 

c) 

d) 


^actual  ”  ^^CR^estimatoc 


h^'^actual 

•  1  ^  ^  10-^*in.-l 

e  ttimated  1 

(17) 

I'^cp'i^  - 

rp%  1 

0.  5% 

6p  %  1 

0.  5% 
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If  any  of  the  conditions,  a  through  d  of  Equation  17,  are  not  satisfied,  steps  Z 
through  5  are  repeated  until  all  errors  are  within  the  specified  limits. 

This  procedure  was  found  to  require  an  average  of  5  trial  computer 
solutions  for  each  operating  point  obtained. 

K.  Estimate  of  Errors 

A  7  X  7  mesh  was  used  in  the  numerical  solution  of  the  Reynolds  and 
Energy  Equations.  This  was  the  finest  mesh  rize  that  could  be  accommodated 
with  an  IBM  650  computer  for  the  present  program.  Previous  experience  of  the 
authors  has  indicated  that  satisfactory  accuracy  can  be  achieved  with  the  7x7 
mesh,  provided  there  are  no  sharp  inflexion  points  in  the  film  thickness  profile. 
As  an  additional  check,  the  calculations  for  one  case  were  repeated  on  a  larger 
computing  machine,  using  a  13  x  13  mesh.  The  results  agreed  with  those 
obtained  using  the  7x7  mesh  within  1%. 

The  error  limits  defined  in  Equation  17  were  set  up  in  order  to  limit  the 
number  of  iterations  required  for  each  solution.  On  the  basis  of  calculations 
carried  out  with  smaller  allowable  errors,  the  effects  of  the  limits  set  in 
Equation  17  are  estimated  to  be: 

Error  in  calculated  maximum  temperature  £  5°F 

Error  in  calculated  minimum  film  thickness  £  0.  0001" 

In  the  calculation  of  the  hydrodynamic  oil  flow  and  the  horsepower  loss, 
additional  errors  are  introduced  in  the  numerical  calculation  of  the  pressure 
gradients  at  the  pad  edges  (see  equations  A-13  and  A- 15).  Particularly  at  high 
loads,  where  the  pad bending^deflections  are  correspondingly  large,  errors  in  the 
calculated  values  of  hydrodynamic  oil  flow  and  horsepower  loss  may  be  as  high 
as  20%. 
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m  RESULTS 

1.  Eight  bMringa  ranging  in  sis*  from  19“  to  SO"  diamoUr  woro  analysed, 
each  at  two  speede  in  die  range  100  to  320  RPM.  These  were; 


BEARINO  SIZE  SPEEDS 

(O.D."  X  LD."1  (RPM) 


19  X  9-1/2 
25  X  12-1/2 
31  X  15-1/2 
37  X  18-1/2 
39  X  19-1/2 
41  X  20-1/2 
45  X  22-1/2 
50  X  25 


160  and  310 
120  and  240 
180  and  320 
180  and  320 
150  and  200 
100  and  200 
100  and  170 
100  and  170 


These  bearings  were  all  geometrically  similar,  with  ths  following  properties: 


L  .  1 
R  *  2 


k  ■  0  85 


0.  154 


•  8p%  =  50 


In  all  cases,  6,  8,  and  10  pad  geometries  were  analysed.  The  results  ars  givsn 
in  Tables  1  through  8  and  plotted  in  ngores  7  tihrou^  46. 

2.  The  ahead  and  astern  bearings  of  Ae  U8S  Barry  (DD933)  and  the  ahead 

bearing  of  the  USS  Norfolk  (DLl )  were  studied  at  tteir  full  speed  condition. 
These  are: 


SIZE 

NUhlBEll 

SPEED 

(aD."xLO.") 

OF  PADS 

(RPM) 

Astern  Bearing  USS  Barry* 

26  X  17-1/2 

8 

160 

Ahead  Bearing  USS  Barry 

31  X  16-1/2 

8 

320 

Ahead  Bearing  USS  Norf^ 

35  X  18-1/2 

8 

170 

The  results  are  given  in  Tables  9  through  11  and  plotted  in  Figures  48  through  56. 


*  In  ths  case  of  the  astern  bearing  of  00933,  Figure  47  was  used  to  determine  Tqj^. 
This  is  because  the  grooves  between  the  pade  of  the  bearing  amounted  to 
approximately  35%  of  the  effective  runner  area,  as  compared  with  about  15%  in  the 
other  bearings.  In  the  absence  of  data  for  this  siss  groove,  Figure  47  was  obuined 
from  Figure  4,  considering  the  groove  temperature  rise  to  be  inversely  proportional 
to  the  extent  of  the  grooves 
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3.  In  order  to  estimete  the  effect*  of  p«d  thickneca,  radial  pivot  location, 
groove  temperature  and  bearing  aiae,  aeveral  additional  calculation* 
were  mad*  varying  theee  parameter*  one  at  a  time.  The  calculations 
were  made  for  the  following: 


Bearing  Sise 

(o.d;  xld;) 

31  X  15-1/2 
31  X  15-1/2 
31  X  15-1/2 
31  X  15-1/2 
31  X  15-1/2 
39  X  19-1/2 
19  X  9-1/2 
75  X  37-1/2 
100  X  50 
19  X  9-1/2 
45  X  27-1/2 
75  X  37-1/2 
100  X  50 


No.  of  Pad* 

6,  8  and  10 
6,  8  and  10 
6,  8  and  10 
6,  8  and  10 
8 
8 

6,  8  and  10 
6,  8  and  10 
6,  8  and  10 
6 
6 
6 
6 


Speed* 

(RPM) 

180  and  320 
180  and  320 
180  and  320 
180  and  320 
180  and  320 
150  and  200 
100 
100 
100 
100 
100 
100 
100 


R 

rp% 

0.130 

50 

0.193 

50 

0.154 

53 

0.154 

47 

0.154 

50 

0.154 

50 

0.154 

50 

0.154 

50 

0.154 

50 

0.130 

50 

0.130 

50 

0.130 

50 

0.130 

50 

^GR 


Per  Figure  4 
Per  Figure  4 
Per  Figure  4 
Per  Figure  4 
130®F 
I30°F 

Per  Figure  4 
Per  Figure  4 
Per  Figure  4 
Per  Figure  4 
Per  Figure  4 
Per  Figure  4 
Per  Figure  4 


The  reeult*  are  given  in  Table*  12  through  21  and  plotted  in  Figure*  57 
through  81. 


4.  A  51  -1/2"  O.D.  X  32"  L  D.  bearing  that  wa*  analyaed  under  aeparate 

contract  with  M.  S.  T.  G.  liG.  Dept,  ha*  al*o  been  included  in  this  report. 
The  geometry  and  operating  condition*  were: 


Bearing  Sise 

Speed* 

(O.  D.  X  L  D. ) 

No.  of  Pad*  (RPM) 

R 

Tqr 

51-1/2  X  32 

10  and  12  200  and  400 

0.  1 17 

52.  6  Per  Figure  4 

The  reeult*  are  given  in  Table  22  and  plotted  in  Figure*  82 
through  86. 
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IV  DESIGN  CHARTS 

In  order  to  facilitate  design  and  selection  of  thrust  bearings,  where  the 
outer  diameter  is  roughly  twice  the  inner  diameter,  the  data  in  Tables  1  through 
H  was  used  to  arrive  at  a  set  of  design  charts.  These  are  given  in  Figures  87 
through  98.  Note  that  in  these  charts  (as  in  the  other  figures  in  this  report) 
solid  lines  represent  data  within  the  range  of  calculations  and  dashed  lines 
indicate  extrapolated  values. 

When  the  oil  film  temperatures  in  Tables  1  through  8  are  plotted,  it  is  seen 

that  both  the  maximum  and  the  average  temperatures,  are,  with  good  accuracy, 

functions  only  of  the  unit  load,  number  of  shoes  and  pad  inlet  temperature.  This 

allows  the  maximum  and  average  temperature  to  be  represented  on  a  single  chart. 

Figure  87.  The  accuracy  of  this  chart,  up  to  T^^^^  ■  2350F,  is  -  5°F.  Above 

T  z  235°F.  the  accuracy  is  -'  10®r. 
max 

The  minimum  film  thickness  is  a  function  of  bearing  size  and  speed  as 
well  as  unit  loading,  number  of  shoes  and  pad  inlet  temperature.  It  is 
represented  here  as  a  function  of  these  variables,  in  the  set  of  nine  charts. 
Figures  88  through  96.  The  accuracy  of  these  charts  is  -  0.0001"  within  the 
calculated  regions.  In  the  extrapolated  regions,  errors  may  be  somewhat  larger. 


The  hydrodynamic  oil  flow  per  pad  is  plotted  in  Figure  97,  as  a  function  of 

V  U  av  \ 

•he  dimensionless  parameter  i  ]  As  was  pointed  out  earlier,  the  oil 


\ 


avg 


B 


flow  calculations  are  subject  to  sigi.xticant  error  (in  some  cases  as  high  as  20%). 
in  part  because  of  the  numerical  approximation  of  the  pressure  gradient  at  the 
pad  edges.  It  is  also  necessary  to  keep  in  mind  that  the  oil  flow  given  by' Figure  97 
is  only  that  which  flows  through  the  clearance  spaces  between  the  runner  and  the 
bearing  pads.  The  total  flow  furnished  to  the  bearing  should  bp  several  tunes 
this  quantity. 


The  friction  horaeoower  loss 


function  of  the  dimensionless  parameteri  avg  avg 


per  Dad  is  plotted  in  Figure  93.  also  as  a 

(/\  «  ^ 


\P 


This  horsepower  loss 


avg 


B 


is  dependent  on  the  calculated  oil  flow  aud  Is  thus  subject  to  lb'*  same  errors.  In 
addition,  it  should  be  noted  that  Figure  98  shows  only  the  horsepower  loss  due  to 
fluid  shear  in  the  oil  film.  There  are  additional  losses  in  the  bearing,  such  as 
those  due  to  turbulence  in  the  oil  grooves. 
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The  example  below  illustrates  the  use  of  the  charts. 


Example:  Compare  the  performance  of  6  8  and  10  pad  geometries  far  a  35"  O.  D. 

X  17-1/2"  I.  D.  bearing  at  a  speed  of  280  RPM  aiid  a  unit  load  c£  650  psi. 

(Oil  2190T,  k  =  0.  85  in  all  cases) 


6  Pad 

8  Pad 

10  Pad 

p  at  h  .  =  0.0010"  (Per  Figures  88,  91  and  94) 

avg  min 

528 

520 

453 

p  at  h  .  =  0.  0008"  (Per  Figures  89,  92  and  95) 

‘^avg  min 

623 

663 

620 

p  at  h  .  =  0.  0006"  (Per  Figures  90,  93  and  96) 

^avg  min  ° 

745 

832 

815 

.  .  h  .  at  p  =  650  psi  (by  interpolation) 

min  avg  i-  '  / 

0. 0075" 

0. 0082" 

0.  0077" 

T  at  650  psi  (Per  Figure  87)  -  °F 

max  ® 

237 

215 

207 

T  at  650  psi  (Per  Figure  87)  ■  °F 

avg 

185 

185 

184 

/^avg  Figure  5) 

1. 8x10-^ 

1.  8x10-^ 

1.  8x10-6 

U  -  2Tr(R-L/2)N  -in/sec 

avg 

385 

385 

385 

B  =  (R  -  L/2)  e„  MR  -  L/2)  2k TT  . 

n 

11.  7 

8.  76 

7.  01 

/  "  avg  Uavg\  ^ 

\Pavg  ®  / 

0.  091 

a  122 

0.  152 

f— — § - ^]xl06  (Per  Figure  97) 

^  ^avg  / 

24.  9 

24.7 

24.  6 

Q  GPM 

0.  98 

0,  73 

a  58 

GPM  (=  n  X  Q) 

'  5. 9 

.  5.  8 

5.  8 

f  X  10^  (Per  Figure  98) 

t 

'  0.  97 

!  1.  01 

1.  12 

f  B  L  p  U 

IIP  - 

6600 

1 

3.  77 

i 

\ 

2.  94 

2.  61 

HP  (=  n  X  HP) 

22.  6 

23.  5 

i  26.  1 

! 
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V  DISCUSSION  AND  CONCLUSIONS 

1.  The  two  principal  criteria  of  thrust  bearing  perforn.a-ice  are  the 

minimum  film  thickness  and  the  maximum  temperature.  The  present 
analysis,  which  was  limited  to  6,  8  and  1 0  pad  bearings,  showed  that: 

a)  For  each  condition  of  operation  (bearing  size,  load  and  speed), 
there  is  an  optimum  number  of  pads,  from  the  standpoint  of  minimum 
film  thickness.  This  can  be  seen  by  comparing  the  design  charts, 
Figures  88  through  96. 

b)  The  maximum  temperature  can  be  decreased  by  increasing  the 
number  of  pads  in  the  bearing.  This  gain  is  greatest  in  the  critical 
high  load  regions  as  shown  in  Figure  87. 

Note  also  from  Figure  87  that  the  maximum  temperature  is  a 
very  sensitive  indicator  of  bearing  load.  This  is  in  contrast  to  the 
oil  temperature  which  is  little  influenced  by  load  changes. 


2.  At  low  loads,  the  minimum  film  thickness  occurs  at  the  inside  radius 
of  the  trailing  edge.  However,  as  the  bearing  load  (and  hence  the  pad 
deformation)  increases,  the  point  of  minimum  film  thickness  moves 
toward  the  pivot,  as  shown  in  Figure  99,  This  figure  shows  that  the 
radial  location  of  the  point  of  minimum  film  thickness  moves  quite 
rapidly  towards  the  center  region  of  the  pad.  It  can  be  concluded 
from  this  that  failures  which  result  from  small  dirt  particles  in  the 
oil  film  are  most  likely  to  occur  near  the  pivot.  This  is  borne  out 
by  experience. 

Figure  99  also  shows  that  the  location  of  the  point  of  minimum  film 
thickness  is  dependent  on  the  pad  subtended  angle.  Thus,  it  moves 
inward  from  the  trailing  edge  most  rapidly  in  the  case  of  the  6  pad 
bearing. 

The  marked  divergence,  at  high  loads,  between  the  minimum  film 
thickness  and  the  film  thickness  at  the  inside  radius  of  the  trailing 
edge  is  also  shown  in  Figure  100.  In  fact  this  figure  shows  that  at 
high  loads,  the  film  thickness  at  the  inside  radius  of  the  trailing 
edge  becomes  almost  insensitive  to  load  changes. 
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The  effect  of  pad  thickness  is  illustrated  in  Figures  101.  Note  that 
there  is  an  optimum  pad  thickness  at  each  specific  load,  from  both  the 
standpoints  of  minimum  film  thickness  and  maximum  temperature. 

At  low  loads,  thinner  pads  are  preferable  for  the  deformation  there 
allows  a  more  favorable  film  shape.  At  high  loads,  on  the  other 
hand,  deformations  become  excessive  and  reduce  load  carrvlns 
capacity,  ’  “ 


Since  propeller  shaft  bearings  are  required  to  operate  under  either 
erection  of  rotation,  the  pivot  location  can  be  varied  only  radially. 
Figure  102  shows  the  effect  of  radial  pivot  location  on  minimum  film 
thickness  and  maximum  temperature,  for  several  values  of  unit 
loading.  Both  these  sets  of  curves  indicate  that  there  is  an  optimum 
pi^vot  location,  that  varies  with  unit  loading.  The  optimum  locations 
obtained  from  the  two  sets  of  curves  are,  however,  different. 

Thus,  from  the  standpoint  of  minimum  film  thickness,  the  optimum 
pivot  location  approaches  the  mean  radius  from  the  outer  clrcum- 
ference,  as  the  unit  loading  increases.  From  the  standpoint  of 
maximum  temperature  on  the  other  hand,  the  optimum  pivot  location 
approaches  the  mean  radius  from  the  inner  circumference,  also 
as  the  unit  load  increases. 


poove  mixing  temperature  plays  a  very  important  role  In 
bearing  performance.  Figure  103  shows  the  reduction  in  load 
^rrying  capacity  that  accompanies  a  rise  in  the  groove  temperature. 
™is  reduction  is  a  major  one,  as  indicated  in  the  following  Uble 
(obtained  from  Figure  103): 


h  ,  " 

min 

T  OF 

GR  ' 

P  psl 

avg 

T 

max 

0.0006 

130 

1030 

220 

0.0006 

158 

800 

235 

Thus,  for  i  constant  minimum  film  thickness,  a  reduction  of  28®F 
in  groove  temperature  achieves  an  increase  of  28%  in  unit  load 
together  with  a  reduction  of  150F  in  the  maximum  temperature! 
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6.  For  a  geometrically  similar  series  of  bearings,  the  unit  loading 

will  increase  with  bearing  size  (at  a  given  angular  speed  and  minimum 
film  thickness)  as  shown,  for  example,  in  Figure  81.  This  is  of 
course  due  to  the  higher  surface  speed  of  the  larger  bearings.  Note 
however,  that  the  slope  of  the  curve  decreases  quite  rapidly  due  to 
the  rise  in  groove  mixing  temperature  and  bending  deflections.  This 
points  up  again  the  importance  of  these  two  factors  on  bearing  per¬ 
formance. 

7.  Early  in  this  report,  it  was  pointed  out  that  the  inclusim  of  thermal 
wedge  and  pad  bending  in  the  analysis  explains  the  load  carrying  capacity 
of  centrally  pivoted  pads.  The  load  carrying  capacities  of  a  flat  pad 
bearing  with  optimum  pivot  location  and  of  a  centrally  pivoted  pad  were 
then  compared  in  Figure  2,  for  a  particular  pad  geometry.  It  should 

be  noted,  howe\er,  that  the  hydrodynamic  pressure  profiles  differ 
markedly  in  the  two  cases,  as  shown  in  Figure  104. 

8.  In  the  present  analysis,  heat  conduction  was  neglected.  Thus,  the 
calculated  maximum  film  temperatures  are  somewhat  higher  than 
those  which  occur  in  practice.  (The  calculated  values  are  therefore 
conservative. )  Furthermore,  whereas  the  calculated  maximum 
temperatures  are  at  the  trailing  edge,  in  practice  they  will  occur  at 
a  small  distance  inward,  also  because  of  conduction. 
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V.  RECOMM  E  NDA  TIONS 

I.  The  ■implified  analyai*  that  was  used  here  has  shown  that  several 

aspects  of  bearing  geometry,  such  as  number  of  pads,  pad  thickness 
and  radial  pivot  location,  have  a  significant  effect  on  load  carrying 
capacity.  The  effect  of  number  of  pads  was  studied  for  a  large  range 
of  bearing  .  izes.  The  effects  of  the  other  factors  were  studied  for 
a  31"  O.  D.  X  15-1/2"  1.  D,  bearing.  It  is  desirable  to; 

a)  Verify  the  results  using  a  more  rigorous  elasticity  analysis  (as 

is  being  done  in  Phase  II). 

b)  Extend  the  results  obtained  to  bearings  with  different  L/R  ratios. 

^  )  Obtain  experimenul  verification  (Phasalll  and  contemplated  thrust 
bearing  test  machine.  ) 

!.  Groove  mixing  temperature  plays  a  very  important  role  in  determining 
the  load  carrying  capacity  of  thrust  bearings.  Gains  in  load  capacity 
on  the  order  of  25%  can  be  achieved  if  mixing  in  the  grooves  can  be 
inhibited,  thus  lowering  the  pad  inlet  temperature  of  the  oil.  This 
suggests  an  investigation  aimed  at  developing  suitable  baffles  in  the  oil 
grooves  which  would  reduce  the  carry  over  of  hot  oil  from  the  downstream 
pad.  (A  reduction  in  the  dl  temperature  at  housing  inlets  also  serves  to 
improve  performance.  ) 

The  major  importance  of  pad  geometry  and  groove  mixing  temperature  on 
bearing  performance  indicate  that  design  modifications  can  be  made  to 
greatlv  increase  the  load  carrying  capacity  of  tilting  pad  bearings.  In 
such  designs,  consideration  should  be  given  to  multi-point  supports  and 
to  shaped  pad  surfaces  as  well  as  to  the  other  aspects  of  bearing  geometry 
studied  in  this  report.  Advantage  should  also  be  taken  of  the  elasticity 
of  the  pads  in  optimizing  the  bearing  design. 

The  present  analysis  was  limited  to  steady  state  conditions.  Analytical 
and  experimental  investigations  are  necessary  in  order  to  arrive  at 
means  of  predicting  bearing  performance  under  transient  conditions, 
such  as  acceleration,  crash  maneuvers,  start  up  under  load  (as  in  a 
submerged  submarine)  and  others. 

Metallurgical  work  is  badly  needed  to-day  to  set  up  operating  temperature 
limits  of  the  various  babbitts  in  use  as  well  as  to  develq>  alternative  materials. 
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6.  In  future  analytical  work,  the  effect  of  thermal  conduction  should  be 
studied. 

7.  The  extent  of  misalignment  present  in  thrust  bearing  installations  on 
board  ship  needs  to  be  investigated.  In  parallel  with  this,  the  degree  of 
load  equalization  between  pads  and  the  load  carrying  capacity  of  the 
bearings  under  misalignment  should  be  analyzed. 


APPENDIX 
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APPEN0IX 


1.  Finite  Difference  Equationa 

The  numerical  solution  of  the  Reynolds  and  Energy  Equations  by  means  of 
finite  differences  is  described  in  Reference  11.  Here  a  brief  outline  of  the 
procedure  is  given. 

For  convenience  of  comparison  between  geometrically  similar  pads,  the 
Reynolds  and  Energy  Equations  are  first  put  in  dimensionless  form. 

Let  r  =  Rf  r2 

h  =  h^h  P  =  12  Inhere  N'»  » 

®  ®  -  121TNV 

"  T  a  -  T 

^g  JCp 

Introducing  Equation  (A-1)  into  Equation  (1)  of  the  text,  we  obtain  the  Reynolds 
Equation  in  dimensionless  form; 


(A- 2) 


(A- 3) 


And,  introducing  Equation  (A-3)  into  Equation  (A-2)  and  solving  for  P.  ,  we  obtain 

*•  J 


i  1 

\  H  1 

1+1/2. J 

H^Cf 

1  Pi  1  ^l.-1+l\ 

li-1/2.1  d04  ) 

(fl 

+  2!l 

i+i/2,j  ]r  1 

(P\  ? 

\iif  li,J-l/2  ■ 

1-1/2,;) 

Z?T  1 

"ij 

(ril|i,j+l/2  fii|i,j-l/2)  aP 

(A-4) 

A? 
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1  p 
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i-l/2,j) 

^  1 

+111  \  1 

Ifp  l,J+l/2  ril|i,j-l/2;  A^ 
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Similarly,  from  Equations  (A-1)  and  Equation  (2)  of  the  text,  we  obtain  the  Energy 
Equation  in  dimensionless  form; 


(A- 5) 


Referring  to  Figure  6,  the  above  equation  can  be  reduced  to  a  difference  equatfon; 


Now  eolving  for  Tj  we  obtain 
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The  pressure  and  temperature  profilea  are  obtained  by  the  numerical  lolution 
of  the  two  seta  of  E4uationa{'A-4)aad^A-7)ttaing  the  iterative  procedure  described 
on  Page  13. 

The  load  carried  by  each  bearing  pad  (Equation  5  of  the  text)  is  obtained 
by  a  numerical  integration  of  the  pressure  field  over  the  pad  surface. 

S  .  IPl.jlk  'i,J  !*-•) 

The  oil  now  out  of  each  pad  is  given  in  Equations  7  of  the  text.  In  order 
to  generalise  the  solutions,  a  dimensionless  factor  is  used; 


Q 


231  Q 

bOtrNRL  h^ 


(A- 9) 


Equations  (a-  1/  andCA-9)are  introduced  into  Equation  7  and  four  point  approximations 
are  used  for  the  pressure  gradients  at  the  pad  edges  (Reference  12).  Noting  that 
the  pressure  zero  along  the  edges,  the  numsrical  form  of  the  flow  equations 
becomes: 


.  m  —  j 


j=i  /" 


5,  =JAr-  <?h., 


(A-10) 


p  ^  (^)i  mtl/2  ■  1-50CP 

L  -^9  j  ,  jjwr  '"♦‘/c  l.m 


V  0.333P  ) 

nr»-l/2  i,m-l 


^5.  =  —  7—  Z  )  ,  (3.000?  .-  1.  500  p  ,-,+  a.33p  ) 

*  L  -^'r  J.i/r  nU/2.  J  "-J  *n-l/2  n-l,j 


The  total  flow  out  of  the  bearing  is  given  by 

=  Q2  +  Q3  +  Q4 


Q 


(A- 11) 
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*•  **®*»*P®wr  lots  by  fluid  ahuar  in  each  baaring  pad  was  written  in 

Eqiation  8  of  the  text  in  the  form 


HP  : 

60  0.  707 


(A- 12) 


In  dimeneionleaa  form  thia  ia: 

HP  =  0  A  r 

where,  from  Equationa  (A-l),  (A-9)  and  <A*12): 
—  .  0.  707 


12  tr^ 

In  numerical  form  the  horaepowar  loaa  la  then  given  by: 
HP,  I  ^ 


(A- 13) 


(A- 14) 


<—  ^  a  _ 
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—  -  n 

R  Ar 
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and  the  toUl  horaepower  loaa  In  each  bearing  pad  ia; 
HP  =  HF2+  H^  +  H^ 
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The  equations  for  the  center  of  pressure  (Equations  10  of  the  text)  are 
numerical  form: 


m 

Ae  Ai- 

_ iii 


n 


z: 


r 

;  , 


sin  6  ^  p  ) 

ij  i.jk 


(A- 16) 
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1.  j  1.  J  ■  J 


W 


n;e; 


Vi  0p% 


100  tan*^  (x/y) 


(A- 17) 
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2.  Film  Shape 


( 

\ 


FIGURE  A-1 


\Ve  consider  that  the  convex  shape  to  which  the  pad  bends  under  load  may  be 
represented  by  part  of  a  spherical  surface  whose  radius  of.curvature  is  R^.,  as  shown 
in  Figure  A-1.  In  all  cases  considered  here,  R^.  is  very  large,  greater  than  10'* 
inches.  In  addition,  the  pad  inclines,  so  that  its  Ungent  plane  directly  above  the 
pivot  point  has  slopes  mg  (circumferentially)  and  m^,  (radially),  with  respect  to 
the  plane  of  the  runner,  as  shown  in  Figure  A-1.  The  pad  inclinations  are  small 
so  that 

sin  mg  -  tan  mg  -  mg 


cos  m  g  =  1 

(A- 18) 

sin  m  =  tan  m_  *  m_ 
r  r  r 


r 


cos  m 


1 
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Let  the  film  thickness  at  a  reference  point  (x^i  y^)  on  the  pad  surface  be  h^. 
The  film  thickness  at  any  other  point  (x,  y)  can  then  he  written; 


h  :  ha  -  mg  (x  -  x^)+  (y  -  y^)  +  R^, 


Since  R|.  is  very  large,  powers  of  the  ratio  /r^  Wxe 
then  becomes:  wj 


2  2\  i 

X  +  y 

-TT- 


(A.19) 


neglected.  Equation  (A- 19) 


h  :  h^  -  mg  (x  -  x^)  +  m^  (y  -  y^)  + 


(x^+  y2)  -  (x2  +  y2) 


2  R- 


(A- 20) 


This  equation  can  be  converted  from  the  x,  y  co-ordinate  system  to  the  r,  6 
coordinate  system  of  Figure  A-1  by  means  of  the  relations: 


X  =  r  sin  (0  -  -i)  -  r_  sin  (0  -  ^  ) 

?  P  P  2 
®  r  ®T 

y  s  r  cos  (0  -  y)  -  Tp  cos  (0^  - 


(A- 21) 


The  general  equation  for  the  film  shape  in  polar  coordinates  is  then: 


h  =  h^+  mg 


a 


sin  (0^  -  -^)  -  r  sin  '0  -  y) 


r  cos  (0 - — ) 

a  a  2 


(  A- 


-  r  cos  (6 - —  ' 


2  R 


r^  -  -  2rrp  cos  (0-  Op)  t  2r^rp  cos  (0^-  0^) 


(note  that  Equation  A- 22  can  also  be  used  to  describe  the  film  shape  for  flat 
pads.  In  such  cases,  R^  is  infinite,  thus  eliminating  the  fourth  term  on  the  right 
hand  side  of  the  equation. ) 
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3.  Bending  Coefficient 

Equation  A- 20  and  Figure  A-1  show  that  (with  the  simplified  elasticity 
approach  used  here),  the  bending  deflection  along  a  point  on  the  pad  surface  is 
proportional  to  the  square  of  its  distance  from  the  pivot,  i.e. 

f=  K(x2  +  y2)  (A-23) 

where  K  »  l/|^2Rc)  (A- 24) 

The  value  of  the  bending  coefficient  K  was  obtained  by  calculating  the 
deflection  at  the  rim  of  an  equivalent  circular  plate  point  supported  at  the  center 
of  its  lower  face  and  carrying  a  conically  distributed  load  on  its  upper  face. 

A  circular  plate  was  used  because  a  closed  solution  for  its  bending  deflections 
is  available  (Reference  11).  A  conical  load  distribution  was  selected  because  the 
ratio  of  peak  to  average  pressure  (3:1)  is  similar  to  that  in  an  actual  bearing 
pad. 


Integrating  Equation  57  of  Reference  11,  for  a  steel  circular  plate  (radius 
"a"  and  thickness  "fnvo")  under  the  loading  and  support  described  above,  the 
deflection  at  the  rim  is  found  to  be: 

S  =  0.75  X  10-8  JiL£_  ,^.25) 

*svg 

From  equations  A-23  and  A- 25,  the  relation  between  the  bending  coefficient  and 
the  pad  load  is: 

K  :  0.  75  X  10-8—11 


This  relation  was  used  in  all  the  Phase  I  solutions. 


(A-26) 
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